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Abstract 

Summation of the perturbation series for the Gell-Mann-Low function /3(g) of (/> 4 theory 
leads to the asymptotics /3(g) = (3oog a at g —> oo, where a ~ 1 for space dimensions 
d = 2,3,4. The natural hypothesis arises, that asymptotic behavior is /3(g) ~ g for 
all d. Consideration of the ’’toy” zero-dimensional model confirms the hypothesis and 
reveals the origin of this result: it is related with a zero of a certain functional integral. 
This mechanism remains valid for arbitrary space dimensionality d. The same result for 
the asymptotics is obtained for explicitly accepted lattice regularization, while the use 
of high-temperature expansions allows to calculate the whole /3-function. As a result, 
the /3-function of four-dimensional </> 4 theory is appeared to be non-alternating and has 
a linear asymptotics at infinity. The analogous situation is valid for QED. According 
to the Bogoliubov and Shirkov classification, it means possibility to construct the con¬ 
tinuous theory with finite interaction at large distances. This conclusion is in visible 
contradiction with the lattice results indicating triviality of </> 4 theory. This contradic¬ 
tion is resolved by a special character of renormalizability in </> 4 theory: to obtain the 
continuous renormalized theory, there is no need to eliminate a lattice from the bare 
theory. In fact, such kind of renormalizability is not accidental and can be understood in 
the framework of Wilson’s many-parameter renormalization group. Application of these 
ideas to QCD shows that Wilson’s theory of confinement is not purely illustrative, but 
has a direct relation to a real situation. As a result, the problem of analytical proof of 
confinement and a mass gap can be considered as solved, at least on the physical level of 
rigor. 
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1. Introduction 


In 1954 Landau, Abrikosov and Khalatnikov |T| derived the famous relation between the bare charge 
go and observable charge g for renormalizable field theories: 


9 = 


9 o 

1 + ^2 go In A /m ’ 


(1) 


where m is the mass of the particle, and A is the momentum cut-off. The constant 02 is positive in (/> 4 
theory and QED, so g tends to zero in the limit A — >- oo for any finite go, he. the ’’zero charge” situation 
takes place. In fact, the proper interpretation of Eq. 1 was given in [L and consists in its inverting, so that 
go is attributed to the length scale A” 1 and is chosen to give a correct value of g: 

1 — /?25lnA/m ^ ^ 


The growth of go with A invalidates Eq. 1 (obtained perturbatively) in the region go ~ 1, and existence of 
’’the Landau pole” in Eq. 2 has no physical sense. 

A little later, Landau and Pomeranchuk [2] put forward arguments on validity of Eq. 1 for arbitrary go- 
They have noticed that the constant limit for the observable charge g can be obtained in the limit go —!► oo 
from the functional integrals of the lattice </> 4 theory, if the quadratic in </> terms are omitted in the action 
On the other hand, the constant limit 1 /(/?2 In A/m) is reached with growth of go already in the 
weak coupling region, which is described by Eq. 1. It looks that neglecting of quadratic terms is possible 
already for g 0 <C 1 , and it is all the more possible for g 0 > 1 : it gives a reason to consider Eq. 1 to be valid 
for arbitrary go. Analogous arguments are possible in the case of QED [2]. These results lead Landau to 
conclusion on fundamental deficiency of the field theoretical description [3}. 

This conclusion was questioned by Bogoliubov and Shirkov [4], who noted that actual behavior of the 
charge g(L) as a function of the length scale L is determined by the Gell-Mann - Low equation 


~d!nZ = ^ = fog 2 + fag 3 + ■ ■ ■ ( 3 ) 

and depends on appearance of the function /3(g). According to the Bogoliubov and Shirkov classification 
0], there are three qualitatively different possibilities (Fig. 1): (i) if /3(g) has a zero at some point g*, 
then the effective coupling g tends to g* at small L; (ii) if /3(g) is non-alternating and has asymptotic 
behavior g a with a < 1, then g(L ) grows to infinity; (iii) if non-alternating /3(g) behaves at infinity as g“ 
with a > 1, then g(L) is divergent at some finite Lq (the real Landau pole arises) and dependence g(L) 
is not defined at smaller distances: the theory is internally inconsistent and a finite interaction at large 
distances is impossible in the continual limit. The latter case corresponds to the ’’zero charge” situation 
in full theory, beyond its perturbative context. Realization of this situation cannot be proved, because a 
behavior of the /3-function is unknown. 0 

In current literature these problems are discussed in relation with the concept of ’’triviality”, introduced 
by Wilson [5j- I n the theory of critical phenomena, a finite interaction is accepted at small length scales 
corresponding to the lattice spacing, while Eq. 3 is integrated in direction of large L. If /3(g) is positive 
beyond the origin, then g —> 0 at large distances and </> 4 theory reduces to the trivial Gaussian model: it 
corresponds to the absence of interaction between large-scale fluctuations of the order parameter. According 
to Wilson’s renormalization group [B] such triviality takes place for Euclidean </> 4 theory in space dimensions 
d > 4. Success of Wilson’s e-expansion |B] is directly related with this triviality: for d = 4 — e, interaction 
between large-scale fluctuations becomes finite but small for e <C 1. 

1 In fact, the limit go —► oo is accompanied by the limit m 2 —> —oo for the bare mass, so the terms and go</> 4 are 

equally significant (Sec. 6). Accepting different laws of growth for — m^, one can obtain different possibilities. On the other 
hand, the strong coupling limit for g can be attained for complex go (Sec. 4) where arguments by Landau and Pomeranchuk 
are not valid in principle. 

2 Equation (1) follows from (3), if only the first term is retained in the right hand side. It cannot be exact due to finiteness 

of /3 3 . 


2 






g(L) 




Figure 1: Three qualitatively different situations according to the Bogoliubov and Shirkov classification. 


In the weak coupling region, the /3-function of four-dimensional </i 4 theory is positive and triviality surely 
exists. In subsequent papers, Wilson set problem more deeply: does triviality for d = 4 exist only for small 
go, or has the global character? Using logic of proof by contradiction, he assumed existence of the boundary 
gf for the domain of attraction of the Gaussian fixed point g = 0 (which is equivalent to alternating behavior 
for /3(g)) and derived the consequences convenient for numerical verification. According to his results [5], 
there are no indications on existence of gf. Historically, it was the first real attempt to investigate the 
strong coupling regime for c/> 4 theory and the first evidence of non-alternating behavior of (3(g)- 

Another definition of triviality was given in the mathematical papers m-m- It corresponds to true 
triviality, i.e. impossibility in principle to construct continuous theory with finite interaction at large 
distances. It is equivalent to internal inconsistency in the Bogoliubov and Shirkov sense, or Landau’s ’’zero 
charge”. It was rigorously proved in m-m that (j) A theory is trivial for d > 4 and nontrivial for d < 4; 
using experience of these proofs, some plausible arguments were given in favor of triviality for d = 4. From 
the physical point of view, the former results are rather evident m ■■ triviality for d > 4 follows from 
nonrenormalizability of </> 4 theory, while nontriviality for d < 4 is a consequence of the nonzero root of 
/3(g), whose existence is easily established for d = 4 — e with e 1. These results do not require any 
study of the strong coupling region, and hence no propositions can be made for the case d = 4, where such 
investigation is obligatory. 

It should be clear that two definitions of triviality are not equivalent. Wilson triviality needs only 
positiveness of the /3-function for g ^ 0, while true triviality demands in addition the corresponding 
asymptotic behavior. Wilson triviality can be considered as firmly established (see a review and numerous 
references in EH), while evidence of true triviality is not extensive and allows different interpretation. The 
most interesting example is shown in Fig.2: it gives dependencies of the renormalized charge against the bare 
one for fixed A/m |llj and looks as numerical confirmation of argumentation by Landau and Pomeranchuk 
(N is proportional to A/m). More close inspection reveals that all results for finite go correspond to the 
parabolic portion of the /3-function^ and do not manifest essential deviations from Eq. 1, while the points 

3 In </> 4 theory, the ’’natural” normalization of charge corresponds to the interaction term written as (167 t 2 /4!)^0 4 . In this 
case, the nearest singularity in the Borel plane (Sec. 2) lies at the unit distance from the origin, and (3(g) is expected to change 
on the scale of the order of unity. In fact, even in the natural normalization the one-loop behavior appears to be somewhat 
dragged-out, and approximately quadratic dependence of (3(g) continues till g ~ 10 (see Fig.4,a below). If the interaction 
term is written as g<p ^/8 or g</> 4 /4!, the boundary between ’’weak coupling” and ’’strong coupling” regions lies at g ~ 10 3 
instead of g ~ 1. 
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Figure 2: The renormalized charge gn( 0) (estimated for zero momenta) against the bare charge go (corre¬ 
sponding to interatomic spacing a) in four-dimensional </> 4 theory for fixed values of Na and to but different 
number TV 4 of lattice sites (according to jT]). 


for go = oo were obtained by reducing to the Ising model, which is an ambiguous procedure. 

In fact, two definitions of triviality were hopelessly mixed in the literature (see Sec. 8). As a result, to 
the end of 20-th century a conviction in triviality of </> 4 theory and QED became predominated in literature. 
Below we overview the comparatively new results m-m obtained after 2000, which prove the absence 
of true triviality for these theories. 

It is clear from preceding discussion that solution of the ’’zero charge” problem needs calculation of 
the Gell-Mann - Low function /3(g) at arbitrary g , and in particular its asymptotic behavior for g —> oo. 
Approaches to this problem are discussed in the next sections. Summation of the perturbation series for 
/ 3(g) with the use of the Lipatov asymptotics gives the positive /9-function in four-dimensional (/> 4 theory 
and its asymptotic behavior g a with a « 1 (Sec. 2). The same result for a is obtained in dimensions 
d = 2 and d = 3. The arising hypothesis (3(g) ~ g for the asymptotic behavior is confirmed in the zero¬ 
dimensional case (Sec. 3) and extended to arbitrary dimensions (Sec. 4). The same approach allows to 
obtain the /3-function in QED (Sec. 5). The problem of complex-valuedness of the bare coupling constant is 
discussed in Sec. 6 and a scheme without complex parameters is formulated. The latter involves the explicit 
lattice regularization and reveals the surprising property in renormalizability of c/> 4 theory: the continual 
limit in the renormalized theory does not demand the continual limit in the bare theory. This property 
allows to give a final solution of the triviality problem (Sec. 8) and makes it possible to use high-temperature 
expansions for calculation of the /3-function with good precision (Sec. 7). The character of renormalizability 
discovered for c/> 4 theory is shown to have a general character (Sec. 9), which can be applied to justification 
of the Wilson theory of confinement (Sec. 10). 

2. Summation of perturbation series 

Let consider the typical problem in field theory applications. A certain quantity W(g) is given by its 
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formal perturbation expansion 

OO 

W{g) = Y W N (-g) N (4) 

N =0 

in powers of the coupling constant g. The coefficients Wn are given numerically and have the factorial 
asymptotics at N —> oo, 

W% s =coa»r(N + b 0 ), (5) 

which is a typical result obtained by the Lipatov method [2D]. We want to find W(g) for arbitrary g while 
the radius of convergence for (4) is zero. 

The standard summation procedure is based on the Borel transformation: each term is divided and 
multiplied by N\, the factorial in the numerator is replaced by the definition of the gamma-function, then 
summation and integration are interchanged, 


OO OO 

OO OO jr-r « p OO 

W(g) = Y W "9 N = E wf / dxx N e~ x g N = / dxe~ x Y ' 

JV=0 N=0 1 ' n n iV=0 


Wn 

m 


(gx) 


N 


and we have a series with a factorially improved convergence. We can use T(N + b) with arbitrary b instead 
TV! and obtain the general Borel-Leroy transformation: 


OO 

W(g) = J dxe~ x x b ~ 1 B(gx) 


B(z) = ZB N (-z)", = ( 6 ) 

The function W{g) is related with its Borel transform B(z) by some integral transformation, while B(z) 
is given by a series with a factorially improved convergence. 

It is easy to show that the Borel transform B{z) has a singularity at the point z = — 1/ao (Fig. 3, a) 
determined by the parameter ao in the Lipatov asymptotics (5). The series for B{z ) is convergent in the 
disk \z\ < 1/ao, while we should know it on the positive semi-axis, in order to perform integration in the 
Borel integral (6); so we need analytical continuation of B(z). Such analytical continuation is easy if the 
coefficients Wjv are defined by some simple formula, but it is a problem when they are given numerically. 

The elegant solution of this problem was given by Le Guillou and Zinn-Justin in 1977 m- It is based 
on the hypothesis that in field theory applications all singularities of B(z) he on the negative semi-axis. 
This hypothesis can be proved in the case of </ 4 theory [22]. □ If such analytical properties are accepted, we 
can make a conformal transformation z = f(u ), mapping the complex plane with the cut (Fig. 3, a) into 
the unit disk |u| < 1 (Fig. 3, b ). If we re-expand B{z) in powers of u, 

OO OO 

B(z) = Y Bn(-z) n \ z = f[u) —► B(u) = Y Unu n , (7) 

N=0 N—0 

then such series will be convergent for any z. Indeed, all singular points P, Q, P,... of B(z) lie on the 
cut, and their images P, Q, Q', P, P',... in the u plane appear on the boundary of the disk \u\ = 1. The 
re-expanded series in (7) is convergent in the disk \u\ < 1, but the interior of the disk is in the one-to-one 
correspondence with the analyticity domain in the z plane. 

4 Validity of this hypothesis is frequently questioned in relation with possible existence of renormalon singularities 1231 . 
Such singularities can be easily obtained by summing some special sequences of diagrams, but their existence was never 
proved, if all diagrams are taken into account m- The given below results for the asymptotics of (3(g) (Secs. 4, 5) are in 
agreement with a general criterion for absence of renormalon singularities m and a proof of their absence for theory m 
(see a detailed discussion in EH). 
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Figure 3: (a) The Borel transform B(z) is analytical in the complex plane with the cut (—oo, —1/ao); 

(6) Its domain of analyticity can be conformally mapped to a unit disk in the u plane; (c) If analytic 
continuation is restricted to the positive semi-axis, then a conformal mapping can be made to any domain, 
for which the point u = 1 is the nearest to the origin of all boundary points; (d) An extreme case of such 
domain is the u plane with the cut (1, oo). 


Such conformal mapping is unique (apart from trivial modifications), if we want to make analytical 
continuation to the whole domain of analyticity. In fact, such strong demand is not necessary since we 
need B{z) only at the positive semi-axis, in order to produce integration in (6). If we accept that the 
image of z = 0 is u = 0 and the image of z = oo is u = 1, then we can make a conformal mapping to any 
domain, for which the point u = 1 is the nearest to the origin of all boundary points (Fig. 3, c). The series 
in u converges for |it| < 1 and particularly in the interval 0 < u < 1, which is the image of the positive 
semi-axis. 

Such kind of conformal mapping has advantage in the strong coupling region. A divergency of the 
series in u is determined by the nearest singular point u = 1, which is an image of infinity: so the large 
N behavior of the expansion coefficients Un is related with the strong coupling asymptotics of W(g). In 
order to diminish influence of other singular points P, Q , Q' ,..., it desirable to move away these points as 
far, as possible. Thereby, we come to an extremal form of such conformal mapping, when it is made on the 
whole complex plane with the cut (l,oo) (Fig. 3,d). Mapping of the initial region (Fig. 3, a) to the region 
of Fig. 3, d is given by a simple rational transformation 

u 

o 0 (l - u) ’ 

for which it is easy to find the relation of Un and Bn, 


£ 

II 

to 

o 

v» = £ 

K—l a ° 

K C%Z\ (N> 1), 

(8) 

where Cft = N\/K\(N — K)\ are the binomial coefficients. 

If W(g) has a power law asymptotics 



W(g) = Wcog* , 

9 -t oo, 

(9) 

then the large order behavior of Un 

U N = UooN “- 1 , 

N —> oo , 

(10) 
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( 11 ) 


00 aQr(a)r(6 + a) 

is determined by the parameters a and Woo- Consequently, we come to a very simple algorithm: the 
coefficients Wn of the initial series (4) define the coefficients Un of re-expanded series (7) according to 
Eqs. 6,8, while the behavior of Un at large N (Eqs. 10,11) is related with the strong coupling asymptotics 
(9) of W(g). 

If information on the initial series (4) is sufficient for establishing its strong coupling behavior (9), then 
summation at arbitrary g presents no problem. The coefficients Un are calculated by Eq. 8 for not very 
large N, and then they are continued according to their asymptotics (10). Consequently, we know all 
coefficients of the convergent series (7) and it can be summed with required accuracy. 

One can apply the described procedure to the perturbation series for the /3-function, 


fag) — fa g 2 + fag 3 + ■ ■ ■ + fag L + ■ ■ ■ + coclqT(n + fa)g N +... , ( 12 ) 

having in mind that several first coefficients (till fa) are known from diagrammatic calculations and their 
large order behavior is given by the Lipatov method. The intermediate coefficients can be found by 
interpolation, the natural way for which is as follows. It can be shown that corrections to the Lipatov 
asymptotics has a form of the regular expansion in l/N: 

(3n = coo,qT(N + fa) jl + + ■ ■ ■ + | • (13) 

One can truncate this series and choose the retained coefficients Ak from correspondence with the first 
coefficients fa,... , fa', then the interpolation curve goes through the several known points and automati¬ 
cally reaches its asymptotics. To variate this procedure, one can re-expand the series (13) in the inverse 
powers of N — N, 


(3n 


c 0 a^T(N + b 0 ) 



M 

N-N 


A-2 

(N - N) 2 


A k 

(N-N) 



(14) 


and obtain a set of interpolations, determined by the arbitrary parameter N. 

In the case of four-dimensional (/> 4 theory, a realization of this program m gives the non-alternating 
/3-function (Fig. 4, a), with the results for the exponent a shown in Fig. 4, b. The exponent a is practically 
independent on N, and only its uncertainty depends on this parameter. If we take the result with the 
minimal uncertainty, we have a value a = 0.96 ± 0.01, surprisingly close to unity. @ 

Something close to unity is obtained also in two and three dimensions [2H 30 (Fig. 5). The natural 
hypothesis arises, that fag) has the linear asymptotics 

fag) ~ 9 > S-too (15) 


for arbitrary space dimension d. If this hypothesis is correct, then there is a natural strategy for its 
justification: 

(i) to test it in a simple case d = 0; 

(ii) to find out the mechanism leading to this asymptotics; 

(iii) to generalize this mechanism for arbitrary d. 

Surprisingly, this program can be realized and Eq. 15 is our main result. Since summation of the series 
gives non-alternating fag) (Fig. 4, a), we may conclude that the second possibility of the Bogoliubov and 
Shirkov classification is realized. 


5 Estimation of errors was made in a framework of a certain procedure worked out in m- Subsequent applications have 
shown that such estimation is not very reliable. 
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Figure 4: (a) General appearance of the /3-function in four-dimensional </> 4 theory according to ['1'2] (solid 
curve), and results obtained by other authors (upper, middle, and lower dashed curves correspond to 
[26l l27l [28] respectively). ( b ) Different estimations of the exponent a according to [T2j . 
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Figure 5: Estimations of the exponent a for c/> 4 theory in two and three dimensions mum. 











3. ”Naive” zero-dimensional limit 


Consider the 0(n)-symmetric </> 4 theory with the action 

= f d d xl± £(V<M 2 + \ m l + fe 

I ct a \ a. 

u = goA e , e = 4 — d (16) 

in d-dimensional space; here too is a bare mass, A is a momentum cut-off, go is a dimensionless bare charge. 
It will be essential for us, that the /3-function can be expressed in terms of the functional integrals^. The 
general functional integral of (/> 4 theory 

Z L A ila M ( Xl ’-- *i x m ) = J D(j) (j) ai {xi)(j) a2 {x 2 ) ■ ■■ (t> aM {x M ) exp (-S {(/)}) (17) 

contains M factors of <j> in the pre-exponential; this fact is indicated by subscript M. 

We can take a zero-dimensional limit, considering the system restricted spatially in all directions. If its 
size is sufficiently small, we can neglect the spatial dependence of (f>(x) and omit the terms with gradients 
in Eq. 17; interpreting the functional integral as a multi-dimensional integral on a lattice, we can take the 
system sufficiently small, so it contains only one lattice site. Consequently, the functional integrals transfer 
to the ordinary integrals: 

= J d n <t>K ! • • • 4>ocm exp - |u</> 4 ) . (18) 



This is the usual understanding of zero-dimensional theory. Such model allows to calculate any quantities 
with zero external momenta. If external momenta are not zero, the model is not complete: it does not 
allow to calculate the momentum dependence. To have a closed model, we can accept that there is no 
momentum dependence at all0. This ’’naive” model is internally consistent but does not correspond to 
the true zero-dimensional limit of </> 4 theory^. The latter fact is not essential for us, since this model is 
used only for illustration and the proper consideration of the general d-dimensional case will be given in 
the next section. 

Expressing /3-function in terms of functional integrals, we obtain it in the form of parametric represen¬ 
tation 


P 


9 = 1- 


n K4K0 


n - 


2 K\ 


2n K±Ko 
i + 2 A'| 


k 6 i< o _ 1 
K 4 K 2 1 


1 - 


K 4 K 0 


(19) 

( 20 ) 


The right hand sides of these formulas contain the integrals 


/*00 

I<M{t) = / (j) M+n - 1 d(f) exp (- tq i 2 - </> 4 ) , 

Jo 



( 21 ) 


obtained from (18) by simple transformations. According to (19,20), the quantities g and /? are functions 
of the single parameter t: excluding t we obtain the dependence /3(g)- 

Investigation of (19,20) for real t shows that g and /? as functions of t have a behavior shown in Fig. 6, a; 
a combination of these results shows that /3(g) behaves as in Fig. 6, b. We see that variation of the parameter 

6 Definition of the /3-function depends on the specific renormalization scheme. We accept renormalization conditions at 
zero momenta (see Sec. VI. A in ED). so the length scale L in Eq. 3 corresponds to m 1 . 

7 This point is essential in a definition of the Z-factor, which should be chosen so as to give a dependence p 2 with the 
unit coefficient in the denominator of the Green function G(p). In the described ” naive” theory we accept Z = 1, since the 
momentum dependence is absent. 

8 It will be clear below from Eq. 49 that ’’naive” theory is correct for d = 0, while the physically interesting limit of small 
d is singular and leads to the qualitatively different results. 
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Figure 6: (a) Dependence of g and /3(g) on the parameter t. (b) A resulting appearance of /3(g). 

t along the real axis determines (3(g) in the finite interval 0 < g < g* where 


To advance into the large g region, we should consider the complex values of t. 

It appears, that in the complex t plane we should be interested in zeroes of the integrals Am(£)- The 
origin of these zeroes is very simple. There are two saddle points in the integral AAf(t) , the trivial and 
nontrivial, 

<^ci = 0, 4>c2 = \/—t/2 , (23) 

and K M (t) can be presented as a sum of two saddle point contributions: 

I< M (t) = Aj* +A 1 e i ' pl . (24) 

If these two contributions compensate each other, then the integral can turn to zero. Such compensation 
can be obtained by adjustment of the complex parameter t = \t\e lx , and in fact there are infinite number 
of zeroes lying close to lines % = ±37r/4 and accumulating at infinity (Fig. 7). The above saddle-point 
considerations can be rigorously justified for zeroes lying in the large |i| region. In fact, it is only essential 
for us that (i) zeroes of AAf(t) exist in principle, and (ii) zeroes of different integrals he in different points. 

Now return to the parametric representation (19, 20). It appears, that large values of g can be achieved 
only near the root of the integral A^. If K 2 tends to zero, then (19,20) are simplified, 


n K 4 Kq 
' n + 2 A'f ’ 


P(g )« — 


4 n K 4 K 0 
i + 2 K$ ’ 


and the parametric representation is resolved in the form 

P(g) = 4 g, 00 . 

We see that, indeed, the asymptotic behavior of (3(g) appears to be linear. 


4. General d-dimensional case 
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Figure 7: Zeroes of the integrals A'm( 0 in the complex t plane. 


The same ideas can be applied to the general d-dimensional case. First of all, the actual functional 
integrals can turn to zero by the same reason. Indeed, the complex values of t with large |t| correspond 
to complex go with small | < 70 1 (see Eq. 21), and we come to miraculous conclusion: large values of the 
renormalized charge g corresponds not to large values of the bare charge go (as naturally to think^]), but 
to its complex values; more than that, it is sufficient to consider the region |c/o| <C 1, where the saddle- 
point approximation is applicable. As a result, the zeroes of the functional integrals can be obtained by 
compensation of the saddle-point contributions of the trivial vacuum and of the instanton configuration 
with the minimal action; contributions of higher instantons are inessential for | < 70 1 "C 1. 

Now we need representation of the /3-function in terms of functional integrals. The Fourier transform 
of (17) will be denoted as Km after extraction of the S -function of momentum conservation and a factor 
Iai...a M depending on tensor indices: 

Z k = K M(Pi) NS pi + . ■■+PM (27) 


where N is the number of sites on the lattice, and I ai ... aM is a sum of terms like 5 aia2 5 a30Ci ... with all 
possible pairings. In general, integrals Km(Pi) are taken at zero momenta, and only the integral K 2 should 
be known for small momentum 


k 2 (p) = k 2 - k 2 p 2 + ... 


(28) 


Expressing the /3-function in terms of functional integrals, we have a parametric representation 
for details): 


(K 2 \ d/2 K 4 Ko 

9 UJ n 


(see US] 


(29) 


(K£\ d/2 / , K 4 Kq o (K' 4 Ko + K a K' 0 )K 2 - 2K 4 KqK' 2 k 2 
\k 2 ) \ Kl K'i K 2 k! 2 - K! 2 k 2 


(30) 


9 It is commonly accepted that the bare charge go is the same quantity as the renormalized charge g at the length scale 
A -1 . In fact, these two quantities coincide only on the two-loop level m and this relation is valid only if g <C 1 and 30 CI 
simultaneously. 
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Figure 8: General appearance of the /3-function in QED 1131 . 


If go and A are fixed, then the right hand sides of these equations are functions of only too, while dependence 
on the specific choice of go and A is absent due to general theorems [31]. 

We see from Eq. 29 that large values of g can be obtained near the root of either K- 2 , or A' 2 . If K 2 —> 0, 
equations (29, 30) are simplified, so g and /3 are given by the same expression apart from the factor d, 


(K 2 \ d/ 2 KiKo 

9 UJ K$ ’ 

and the parametric representation is resolved as 

P(g) = dg, 


p = -d 


(K 2 \ d/2 K 4 K 0 

UJ Kl 


g -> 00. 


(31) 


(32) 


For K 2 —>• 0, the limit g —>• 00 can be achieved only for d < 4: 

P{g) = (d-4)g, g^oo. (33) 

The results (32), (33) correspond to different branches of the analytical function /3(g). It is easy to 
understand that the physical branch is the first of them. Indeed, it is well known from the phase transitions 
theory that properties of </> 4 theory change smoothly as a function of space dimension, and results for 
d = 2,3 can be obtained by analytic continuation from d = 4 — e. According to all available information, 
the four-dimensional /3-function is positive, and thus has a positive asymptotics; by continuity, the positive 
asymptotics is expected for d < 4. The result (32) does obey these demands, while the branch (33) does 
not exist for d = 4 at all. Eq. 32 agrees with the approximate results discussed in Sec. 2 and with the exact 
asymptotic result /3(g) = 2 g, obtained for the 2D Ising model [32, from the duality relational. 


5. Calculation of /3-function in QED 

The same ideas can be applied to quantum electrodynamics. Summation of the perturbation series for 
QED 13) gives the non-alternating /3-function (Fig. 8) with the asymptotics /3oog a , where (Fig. 9) 

a =1.0 ±0.1, /3oo = 1.0 ± 0.3 (34) 

(g = e 2 is the running fine structure constant). Within uncertainty, the obtained /3-function satisfies 
inequality 

0 <P(g)<g, (35) 

10 Definition of the /3-function in E2 differs by the sign from the present paper. 
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Figure 9: Different estimations of the parameters a and (3^ for QED according to US- 


established in [341135] from the spectral representations, while asymptotics (34) corresponds to the upper 
bound of (35). Such coincidence is hardly incident and indicates that asymptotics /3(g) = g is an exact 
result. We shall see below that it is so indeed. 

The general functional integral of QED contains M photon and 2TV fermionic fields in the pre-exponential, 

Im,2n = J DADipDil;A IJ , 1 (xi)...A IJ , M (xM)tp(yi)'^(zi)..-ip(yN)i’(zN)exp(-S{A,iJ},'4i}) , (36) 


where S{A,'ip,ip} is the Euclidean action, 


S{A,i(),il>} 


d A x 



V 


dvA^) 2 + ip(i$ — mo + , 


(37) 


while eo and mo are the bare charge and mass. The Fourier transforms of the integrals Im,n with excluded 
(5-functions of the momentum conservation will be referred as KMN(li,Pi ) after extraction of the usual 
factors depending on tensor indices^; g,; and pi are momenta of photons and electrons. 

In general, these functional integrals are taken for zero momenta, but two integrals Kq 2 (p) and K 2 o(q) 
should be estimated with the lowest order momentum corrections: the first is linear in p , and the second 
is quadratic in g, 

Ko2(p) = K 02 + Ko2jf > K2o(q) = K 20 + K2oq 2 , (38) 


and in fact the tilde denotes their momentum derivatives. 

Expressing the /3-function in terms of functional integrals (see [16] for details), we have a parametric 
representation 


K 2 r 2 Koo 
K 2 m K 20 ’ 


(39) 


A specific form of these factors is inessential, since the results are independent on the absolute normalization of e and m. 
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( 40 ) 


= 1 K 02 K 02 K 2 2 K 00 f 2K± + Ko_ 2K^ _K *o 1 

[9> 2 A- 02 /?' 2 - K' 02 K 02 K 2 2 K 20 \ K 12 Koo K 02 K 20 J 

According to Secs. 3,4, the strong coupling regime for renormalized interaction is related with a zero of a 
certain functional integral. It is clear from (39) that the limit g — > oo can be realized by two ways: tending 
to zero either K 02 or K 20 . For K 02 —>• 0, equations (39,40) are simplified, 


KhKoo 

K% 2 K 20 ’ 


0(s) 


and the parametric representation is resolved in the form 


K^Kpo 

Ki 2 K 20 ’ 


(41) 


0(g) =9, 5 -too. 


(42) 


For K 2 o —> 0, one has 


P{g)ocg 2 , j-> oo. 


(43) 


Consequently, there are two possibilities for the asymptotics of /3(g) , either (42) or (43). The second 
possibility is in conflict with inequality (35), while the first one is in excellent agreement with results (34) 
obtained by summation of the perturbation series. In our opinion, it is a sufficient reason to consider (42) 
as an exact result for asymptotical behavior of the /3-function. It means that the fine structure constant 
in pure QED behaves as g ot L~ 2 at small distances. 


6. A scheme without complex parameters 

Our use of the complex bare parameters may look suspicious, since it corresponds to the non-Hermitian 
bare Hamiltonian; at first glance, it violates unitarity, since the S’-matrix is expressed through the Dyson 
T-exponential of the bare action. 

In fact, a problem is solved by Bogoliubov’s construction of the axiomatical S-matrix [J: according 
to it, the general form of the A-matrix is given by the T-exponential of iA , where A is a sum of (i) the 
bare action, and (ii) a sequence of arbitrary ’’integration constants” which are determined by quasi-local 
operators. In the regularized theory we can set the ’’integration constants” to be zero, and the A-matrix 
is determined by the bare action. However, in the course of renormalization these constants are taken 
non-zero, in order to remove divergences. These non-zero ’’integration constants” can be absorbed by the 
action due to the change of its parameters. As a result, for the true continual theory the A-matrix is 
determined by the renormalized Lagrangian, which is Hermitian for real g. 

Nevertheless, scientific community has a bias against complex bare parameters: it is related with the 
old discussion between Lee and Pauli, moderated by Heisenberg, on the exactly solvable model suggested 
by Lee [36). After paper f37], the Lee model was considered as unsatisfactory due to existence of ” ghost 
states”, and this point of view was included in many textbooks. Quite recently [38] it was found that 
this point of view is incorrect and the Lee model is completely acceptable. A key idea of [35] is that 
the complex-valued Hamiltonian can be made Herminian by modification of the inner product for the 
corresponding Hilbert space. 

Below we can calm sceptic spirits and suggest a scheme without complex parameters. In this case we 
accept explicitly the lattice regularization and take the action in the form 

S{<i>} = \a d Y ^x-x'0x^x' + \m 2 0 a d Y& + isoa 2d-4 Y & > ( 44 ) 

x,x' X X 

where we accept A = a -1 (a is a lattice spacing) and restrict ourselves by the case n = 1. Making a change 
of variables 

0 t </> (g 0 a 2d ~ 4 / 4) 1/4 (45) 
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and setting t = (1/go) 1 / 2 as before, we can write the functional integral (17) in the form 


Z (M) {xJ = (2i) 


a r+M 



exp 


,/ x 


t / — trriQ i 


-£« 


(46) 


We accept a = 1 measuring J x -x' and in units of A 2 . We consider t as a running parameter of the 
parametric representation and investigate a singularity at t —> 0, which has a simple origin. For go ^ 1, 
Eq. 46 allows expansion over the gradient term t J x _ x '^ x ^ x '. In zero order in t the integral Z^ has a 
d-functional form in the coordinate representation, x') ~ (5 XX /, and its Fourier transform has no 

momentum dependence; the latter appears only in the first order in t. As a result, the integral K 2 in 
expansion (28) is small in comparison with K 2 , i.e. K 2 /K 2 ~ 1/t, which leads to a singularity at t = 0 in 
(29, 30). This singularity is more complicated than other singularities in the t plane (Fig.7) and needs the 
accurate investigation. 

Expansion of the functional integral (46) in powers of t expresses it in terms of ordinary integrals 


/ OO 

d4> 4> 2k exp {-tml4> 2 - (/> 4 } 

-OO 


(47) 


and parametric representation (29,30) reduces to the form 


n J 0 
9 ~ 1 2t I 2 


d/2 


hh 


1 - 


I4I0 


i+2 I 2 


/3(g) 


= d + 2 


6-*2 

~T r 

1 o 


2 Ij 

T 2 


2 

I 2 h 

I 3 


h 

lo 


i+2I I 

n 1$ 


T 2 

1 2 

r 2 

2 0 


h 

lo 


(48) 


This result can be written in a simple form, if the functions g Z ero(t) and f3 zero (t) are introduced, which 
correspond to the zero-dimensional case and have appearance shown in Fig.6,a: 


9 = 


n_ h 
21 I 2 


d/2 


9zero(j'^ , o) i 



Resolving the parametric representation in the limit t — > 0, we come to the asymptotics 


/3(g) 


d + 


>( 0 ) 


9zero(0)_ 


9 ~+ 00 


(49) 


reducing to the result /3(g) = 2.29 g of the paper 39} after substitution of numerical values. However, 
this result is not final. Instead of the limit t —> 0 for fixed mo one can consider a limiting transition 
under condition fmg —> const with different values of const ; then the general structure of theory remains 
unchanged, but asymptotical behavior of /3(g) will be different. The question arises on the correct character 
of the limiting transition, corresponding to the strong coupling regime. 

The gradient expansions of the considered type was exploited in a number of works [39]- [42] . and 
ambiguity of a strong coupling limit was finally realized by their authors. Nevertheless, the correct character 
of the limiting transition was not established till the paper HZ]. In the framework of the parametric 
representation (29,30) the indicated problem accepts the different form. Deficiency of the result (49) 
consists in the presence of two independent parameters t and fm.g. If one of them is excluded in favor of g, 
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Figure 10: The renormalized mass as a function of the bare mass in the strong coupling region. 


then the /3-function depends not only on g but also on m q, while the latter dependence should be absent 
according to general theorems [31]. The question arises on resolving of this contradiction. Of course, there 
is no real contradiction, because the general theorems suggest that the continual limit A —> oo is already 
taken. Physically it means a fulfilment of the condition 

to 2 <C A 2 , (50) 


which is equivalent to the condition £ a for the correlation length £; it means that a characteristic scale 
of the field variation contains many lattice sites, so further diminishing of a is of no significance. As a 
result of the gradient expansions, the renormalized mass is represented in the form (Fig. 10) 


I< 2 _ n Iq 
~ 


m 2_jA£ = _ - = / 

2 1 I2 


n 0 > 

1 A, 


trriQ » 1 


\tml\ < 1 


(51) 


l/f 2 |TOp| , -imp S> 1, 


and the condition m 2 -C 1 (corresponding to (50) in dimensional units) is satisfied in the case 


tm 0 = —k , 


k > 1. 


(52) 


After replacement c/) 2 —> K(?!> 2 /2 the exponential in (46) accepts a form 


exp 


2 ^^ ^ ^ J'x.—y, 


TI ex p{^ 2 ( 2 44 - 44)} 


(53) 


where the last factor is localized near </£ = 1 and can be replaced by A5(q i 2 — 1). The constant A is 
inessential for a ratio of two integrals and one can set A = 1. As a result, Eq. 46 accepts a form 


Zm{^} = ( tn) 


Jf+M 


n 


. <j) XM exp < -\tn ^2 ■/* 


H^x-1) (54) 


and a functional integral transforms into the Ising sum over values </> x = ±1. In the n-component case one 
obtains a er-model m instead of the Ising model. 
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Now all functional integrals depend on the single variable tK and the right-hand sides of (29,30) are 
functions of only this argument; it defines the /3-function depending only on g. In the physical motivation 
of the limiting transition we implied the conditions providing (50) 

f«l, k»1, tK » 1, (55) 

but in fact only first two inequalities were used in derivation of (54). Therefore, transformation to the Ising 
model is valid under conditions 


t <C 1, k>1, tK is arbitrary. (56) 

In particular, it is valid in the region te <C 1, where gradient expansions are possible and large values of 
the renormalized charge g are reached. Using (51,52), we come to conclusion that strong coupling regime 
of <f > 4 theory corresponds to the limit 

t —> 0 , tm,Q —> —oo , tm 2 —»• 0 , to 2 —> oo , 

so neither to o = const nor m = const is the correct condition for a limiting transition. It means that 
dependencies of Fig.2 are not actual from the very beginning. 

Returning to (49) and setting —tm% = k> 1, we have 

and so the asymptotics of the /3-function is obtained 

P(g) =dg {g ->• oo), (58) 

coinciding with (32). We see that singularity at t = 0 leads to the same result as singularities at the 
complex values of t. 

Representation (54) for functional integrals can be used for calculation of the observable quantities. 
The latter are obtained in the form 

A oba =A dA f A {tK), (59) 

where d A is a physical dimension of the quantity A 0 bs . Using analogous expressions for g and to, 

g = fg(tK), TO 2 = A 2 f m (tK) , (60) 

we can rewrite (59) in the form 

A 0 bs = m dA F(g ), (61) 

which does not contained the bare parameters go, too, A; so (61) gives a ’’theorem of renormalizability” for 
a strong coupling region. It should be stressed that we do not take the continual limit in the bare theory 
and retain the lattice as a convenient instrument for representation of functional integrals, and only the 
lattice spacing a is excluded from the physical results. 

7. Application of high temperature expansions 

Rewriting (56) in dimensional quantities and changing tn —> k, we see that reducing of </> 4 theory to 
the Ising model is possible under conditions 

9o > 1, -gQ 1/2 m 2 0 / A 2 > 1, -g^ml/A 2 = k , (62) 

where k is an arbitrary parameter having a sense of the inverse temperature in the Ising model. Corre¬ 
spondingly, Eqs. 48 have a structure 

g = F(k) , p(g) = Fi(k) , (63) 
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and define the /3-function in the parametric form. At first glance, the condition go /$> 1 corresponds to the 
strong coupling regime and parametric representation (63) is limited only by this regime. However, there 
is another view on this situation. Let strengthen conditions (62) by taking the limit 

go -t oo, -So 1 / 2 m o/A 2 -t oo, -g^ml/A 2 = k (64) 

In this case, transition from (46) to (54) is valid without any approximations and conserves a strict equiva¬ 
lence with the initial c /> 4 theory under a certain choice of its bare parameters. The last property conserves a 
form of the Lagrangian under renormalizations. The taken limit go —)> oo does not mean the same limit for 
the renormalized charge g ; in fact, according to gradient expansions, g varies from infinity to the order of 
unity when k changes from zero to a finite value. Since parametric representation (63) is exact and specifies 
the /3 function in the interval 1 <g < oo, it can be analytically continued and treated as a definition of /3(g) 
at arbitrary values of g. However, a certain doubt arises: does this definition provide the correct results in 
the weak-coupling region? 

An answer to this question can be obtained using high-temperature series (44j . Such series are tradi¬ 
tionally constructed for quantities (superscript c marks the connected diagrams) 

X 2 = ^(0x^o) C , ^2 = x 2 (^x<M C > X 4 = ^2 (^x^y^z^>o) C , (65) 

X X x,y,z 


which coincides up to factors with the ratios K 2 /K 0 , K 2 /K 0 , and K^/Kq of the functional integrals 
specified above; more precisely, 


^ = 2d— = -/ 0 (k) 
K 2 d2 « 


K 2 
I<o 


= 2kX 2 = k/2(«) ■ 


( 66 ) 


K 4 K q _ IX 4 
K 2 ~ 3 xf 


-/4OO, 


where the introduced functions /^(k) will be used below. It was taken into account that there is no zeroth 
term in the expansion of /r 2 in k (see Eq. 68 below), so that all functions /o(k), / 2 (k), and / 4 (k) are regular 
and their expansions begin with the zeroth term. The substitution of ( 66 ) into (29,30) gives 


9 = 



d/2 

Ufa) , 


M = d - [i°/<wr 

g l-«[ln /o(k)]' 


(67) 


The initial parametric representation (29,30) for the /3-function is ’’dead”, since evaluation of functional 
integrals looks hopeless. It becomes ’’alive” in the form (67), since functions fi{n) can be calculated using 
high-temperature expansions!^!. 

For a simple hypercubic lattice with the interaction between the nearest neighbors, the first terms of 
the expansion for functions (65) have the form (in the case d = 4, n = 1) [45] 

X 2 = 1 T 16k: -f- 224 k 2 -(-... 


= 16k + 512k 2 + 33920/3k 3 + ... (68) 

X 4 = -2 - 128k - 4672k 2 

Taking limit k —>• 0, it is easy to obtain the strong coupling behavior for the /3-function and anomalous 
dimensions [18], and then develop their expansion in powers of g ~ 2 ^ d . Below, 14 terms of expansion (68) 
are used, which are given for n = 1 in tables 5, 8, and 11 of the paper | To ■ 

12 The function / 2 (ft) does not enter to (67), but it is actual for calculation of anomalous dimensions m 


18 







The use of Pade-approximants allows to obtain the ,3-function and anomalous dimensions at arbitrary 
g. The general strategy consists in the following. The Ising model has a phase transition in a certain point 
k c = 1 /T c , and a typical physical quantity F(k) has the critical behavior of the form 

Foe (T-T c )~ x oc {k c -k)~ x . (69) 


If F(k) is expanded in k, then a convergence radius of the expansion is limited by the quantity k c ; in actual 
cases, k c is the nearest singularity to the coordinate origin. If the logarithmic derivative of F is taken, 

F' -A 

(lnF)' = — =-1- less singular terms , (70) 

r n — k c 


then the main singularity for it is a simple pole with a residue —A and can be investigated using the Pade 
approximation 0fi] . The Pade-approximant [M/N] is defined as the ratio of two polynomials of degrees M 
and N, 

n n , Pm(k) p 0 +p 1 n + ...+p M K M rrn 

[ ’ Qjv(«) 1 + gi K + ... + g N K N ’ [ ’ 

whose coefficients are chosen to reproduce the first M + N + 1 coefficients in the expansion of (In F)' over 
k. It is known that Pade-approximants successfully predict the nearest singularities of the corresponding 
function if these singularities are the simple poles [441 146] F^l . If A and k c are predicted reliably, then the 
whole function F(re) can be found in the interval 0 < k < k c with a good precision. If such results for /)(«) 
are substituted into the right-hand sides of (67), then the /3-function can be determined in the interval 
g* < g < oo, where g* is the fixed point of the renormalization group. In the four-dimensional case, one 
has g* = 0 and /3(g) is completely determined by the described procedure. 

The use of this strategy in the four-dimensional case is complicated by the existence of logarithmic 
corrections to scaling 03 IS]: 

X 2 ~ r -1 | 1 iit| p , C 2 ~^~T- 1 |lnr| p , X 4 ~ r" 4 | lnr^" 1 , (72) 

a 2 


where r ~ (re c — k) is a distance to the transition and p = (n + 2)/(n + 8). Substitution to (66,67) gives 
for behavior of charge g 

9 = y-j— -r j c 0 = 2/^2 (j 0), (73) 

|1ut| 

where the coefficient of the logarithmic factor is universal. If equations (72) and (73) are fulfilled, then the 
parametric representation (67) automatically reproduces the one-loop result of weak-coupling expansion 
for the /3-function. 

The objective test of Eqs. 72 for lattice models were performed in many works [H]-[:55]. In particular, it 
was convincingly shown in [481 1491 that high-temperature series for the Ising model allow reliable prediction 
of the exponent p. Eq. 73 was confirmed with a satisfactory accuracy in 03 ED- Already these results 
provide the positive answer to the question formulated above: parametric representation (63) gives correct 
results for the /3-function in the weak-coupling region. 

The modified treatment procedure suggested in [18] allows to improve estimation of the constant Co, 
which was not very satisfactory in preceding papers. The relation r = Aq(1 — k/k c ) contains the non- 
universal coefficient Aq which is essential in logarithmic factors. If such factors are extracted from functions 
(65) in accordance with (72) and the Pade analysis is applied to remaining functions, then strong coupling 
estimations of the constant co (Fig.11,a) are close to the theoretical value c Q h = 35.09 but systematically 
exceed it. If the logarithmic factors are extracted in accordance with the next-to-leading logarithmic 
approximation, then estimations of cq improve (Fig. 11, b) and become close to the theoretical value in the 
same range 0.85-^1.06 of the non-universal parameter g (analogous to Ao) where the power law singularities 

13 Usually, one uses diagonal ( M = N) or quasi-diagonal ( M ai N) approximants, whose convergence to the corresponding 
function is proved under the most general assumptions m- 
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Figure 11: The constant cq in Eq. 73 versus Aq in the leading logarithmic approximation (a), and versus g 
in the next-to-leading logarithmic approximation (b). 


of remaining factors are close to theoretical expectations. The exact Co value is realized at g ~ 1.02, and 
this value of g can be used in the further analysis. Extracting from the quantities X2, ^2, X4 all logarithmic 
and power-law singularities and applying the Pade-approximation to remaining regular functions, one can 
calculate the right hand sides of (67) for the whole interval 0 < k < k c and obtain the /3-function at 
arbitrary g (see a solid line in Fig.12). To give an impression on accuracy of the calculation, we show by 
the dotted line the analogous dependence obtained under the assumption of the constancy of indicated 
regular functions, when any information on them is dropped from results. In fact, these regular functions 
are known with the accuracy of several percents. 

Analogously, one can obtain a behavior of the renormalized charge g and the renormalized mass m as 
functions of k (Fig.13), which is essential for characterization of the lattice </> 4 theory. 

8. Is (f> 4 theory trivial? 

In the preceding sections we have established, that the Gell-Mann - Low function /3(g) in four¬ 
dimensional <p 4 theory is non-alternating and has asymptotic behavior /3(g) = 4g at g —> 00. According to 
the Bogoliubov and Shirkov classification (Sec. 1), it means a possibility to construct the continuous theory 
with finite interaction at large distances. This conclusion is in visible contradiction with lattice results 
indicating triviality of (f> 4 theory. 

As we stressed in Sec. 1, one should differ two definitions of triviality. Wilson triviality means that 
integration of equation (3) in the direction of large distances L gives the effective charge g tending to zero; 
this definition implies the massless theory, since in the opposite case the distance scale is saturated by the 
inverse mass. In definition of true triviality one considers the massive theory and suggests finite interaction 
goo for L > to -1 ; a theory is trivial, if integration of (3) in direction of small L gives a divergency at 
finite Lq and does not allow to reach the L —> 0 limit. Such situation is internally inconsistent and means 
incorrectness of the initial suggestion on finite interaction at large distances; in fact, L 0 —> 0 if g^ —> 0. 
Wilson triviality means that /3-function is non-negative and has a zero only for g = 0. True triviality needs 
in addition its sufficiently quick growth at infinity, /3(g) ~ g a with a > 1. According to preceding sections, 
</> 4 theory and QED are trivial in the Wilson sense, but do not possess true triviality. 

Two definitions of triviality were hopelessly mixed in literature m, and there are two main reasons 
for that. Firstly, it is rather difficult to test true triviality in the lattice approach^ 4 !. Secondly, one can 

14 A definition of true triviality in the lattice approach was given in mathematical papers 00. When the lattice spacing 
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Figure 12: A solid line is a calculated function (3(g)/g. The dashed lines are the strong- and weak-coupling 
asymptotic behaviors. A dotted line shows the same results obtained under the assumption of the constancy 
of regular functions. 
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advance arguments that give an illusion of equivalence of two definitions. As illustration to the latter 
point, consider the following reasoning. The only alternative to perturbative approach is to express all 
quantities in terms of the functional integrals. The latter depend on the bare charge go, bare mass mo and 
the ultraviolet cut-off A. Taking into account their dimensional character, one has the following relations 
for the renormalized charge g, renormalized mass m and observable quantities A, 


9 = Fg {go,m 0 /A) , 


m = A F m (g 0 , m 0 /A) , (74) 

A, ; = A di Fi (go, m 0 /A) , 

where d* is a physical dimensionality of A,. Excluding go and mo/A in favor of g and m/A, one has 

Ai = m di Fi (g,m/A) . (75) 

To eliminate the dependence on A we should take the limit m/A —» 0. In the lattice approach, this 
limit corresponds to £/a —> oo, i.e. to the phase transition point. The latter is determined by a zero of 
/3-function, which gives g = 0 in four-dimensional </> 4 theory. 

In this argumentation, Wilson triviality is considered as given, while true triviality is ’’derived" from it. 
Of course, such ’’proof’ is incorrect, because two definitions are surely not equivalent. This shortcoming 
originates from our assumption that a generic situation takes place in Eq. 75, and a limit m/A —> 0 is 
necessary. However, the special case is possible, when the m/A dependence is absent in (75), and a 
limiting transition is irrelevant. This special case fills the "gap” between two definitions and makes them 
not equivalent. 

According to Sec. 6, such special case is actually realized in /> 4 theory. Let return to Eqs. 74 and impose 
the condition m <C A, corresponding to the continuum limit of renormalized theory. If this condition is 
imposed in the region go 1, then </> 4 theory reduces to the Ising model, containing the single parameter 
k, which plays the role of inverse temperature; relations (74) accept the form 

g = F g (k) , 

m = AF m (k) , (76) 

Ai = A di F t (k) . 

So far there is nothing unusual: the condition m/A —> 0 gives a relation between go and mo/A, so all 
functions in Eq. 74 depend on the single parameter, which we denoted as k. The non-trivial point consists 
in the fact that condition m/A <C 1 is sufficient for transformation to the Ising model, but not necessary 
for it. This transformation is possible under the weaker conditions, which are compatible with an arbitrary 
value of m/A (Sec. 6). Excluding k from (76), one obtains the equations 

Ai = m di Fi (g) , (77) 

which are analogous to (75), but do not contain the parameter m/A. As a result, the program of renor¬ 
malization is completely fulfilled, and no additional limiting transitions are necessary. It means that (a) 
we can retain the lattice in the bare theory (as a convenient tool for representation of functional integrals), 
and (b) relation between m and A (or £ and a) can be arbitrary, so a finite value of g becomes possible. 

Usually, the lattice theory contains more parameters than the initial field theory. For example, in 
discretization of the gradient term of (/> 4 theory we obtain a set of the overlap integrals J x , which can be 
chosen rather arbitrary. Then an interesting question arises: if we can retain a lattice in the bare theory, 
then what lattice model should be chosen? 

a tends to zero, the bare parameters go and mo should be considered as functions of a. A theory is non-trivial, if there exists 
some choice of functions go(a) and mo (ft), providing finite interaction at large L; if such functions do not exist, then a theory 
is trivial. Of course, it is rather difficult to test ”existence” or ’’non-existence” in numerical simulations. 
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The answer can be found from Eq. 75. Since dependence on m/A is absent, we can tend this ratio to 
zero. But in this limit (when £/a —>• oo) there are physical grounds for independence of functions F, on 
the way of cut-off. If such independence takes place for m/A —>• 0, it retains for arbitrary m/A due to 
independence of functions Fi on this parameter. In fact, this argumentation implies renormalizability of 
theory (due to which the dependence on A can be excluded) and belonging of the lattice model to the 
proper universality class (inside of which the dependence on the way of cut-off is absent). 

The lattice theory is frequently considered as a reasonable approximation to the true field theory. In 
this case we should accept the condition £ a, which signifies that one has a lot of lattice sites on the 
characteristic scale of variation of field. This condition can be strengthen till £/a —> oo or liberalized till 
£ > a. The first case corresponds to the point of phase transition and gives g = 0. In the second case we 
obtain restriction g < 1 (for the proper charge normalization |18j j. which can be used to obtain the upper 
bound on the Higgs mass [5SJ 1HS] ■ 

In fact, the lattice theory should not be considered as any approximation to field theory, though it is 
possible for go <C 1. The true field theory is continuous from the very beginning and does not contain any 
lattice. The lattice is present only in the bare theory, which is an auxiliary construction and is completely 
removed later. No physical requirements, like £ ~3> a, are relevant for it. If one removes the condition 
£ a, then any values of g become admissible^. In fact, a real designation of the bare theory is to 
represent the relations between physical quantities in the parametric form (74). Such representation has 
no deep sense already due to its ambiguity: it can be written in many different forms, changing go and 
mo/A by any other pair of variables. 

We see that contradiction between the continual and lattice approaches is resolved by a special character 
of renormalizability in c/ 4 theory: correct relations (77) between physical quantities can be obtained for 
the arbitrary value of the parameter a/£, while a dependence on this parameter is absent; to obtain the 
continuous renormalized theory, there is no need to eliminate a lattice from the bare theory. 


9. General situation in renormalizable theories 

The interesting question arises: is such kind of renormalizability related with the specific properties of 
</> 4 theory, or it is a manifestation of some general mechanism? 

We shall see below that the second variant is correct. It can be understood in the framework of 
Wilson’s many-parameter renormalization group (RG) [5]. According to it, the parameters pi of some 
lattice Hamiltonian are considered as functions of the length scale Z.l 16 l The flow of these parameters is 
determined by the RG equations, which can be written in the differential form 


dpi 

d\n{l/a) 


Fi{pk} ■ 


These equations can be linearized near the fixed point 


Pi{l)=p* (for all l) 


(78) 


(79) 


and investigated by the standard methods of linear algebra. The ordinary phase transitions are described 
by the saddle points of such equations. The simplest saddle point in the two-parameter space has the 
straight-line trajectories in two main directions (one stable and one unstable), while the rest of trajectories 

15 This point of view is in complete agreement with mathematical definitions El El, according to which the limit a — > 
0 is taken for the arbitrarily chosen dependencies go(a) and mo (a) (see Footnote 14). We impose conditions go —»• oo, 
g Q 1 ^ 2 m 2 a 2 — >• —oo, gQ 1 m 2 a 2 = — k, necessary for transformation to the Ising model (Sec. 6). 

16 Physically it is explained by the well-known Kadanoff construction. In the description of magnetics, one begins with 
the microscopic Hamiltonian for elementary spins in the lattice sites. Then it is possible to introduce the macroscopic spin 
variables corresponding to the blocks of size l and write the effective exchange Hamiltonian for them. Since the blocks of size 
nl can be composed of n d blocks of size l, then recalculation Pi(l) Pi(nl) is possible, i.e. pi(nl) = Hi (n,{pfc(Z)}). Taking 
n close to unity, one can obtain Eqs. 78. 
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Figure 14: Schematic of the Wilson many-parameter space. 


are hyperbolic. For the usual phase transitions, there are infinite number of stable directions and one (in the 
simplest case) unstable direction. The latter is related with some controlling parameter like temperature, 
measuring the distance to the critical point. 

Instead of increasing l for fixed a, we can diminish a for fixed l. The continuum limit a — > 0 of field 
theory corresponds to the critical surface £/a = oo in the many-parameter space (Fig.14). All trajectories 
on the critical surface tend to the fixed point. The unstable trajectory, originating at the fixed point will 
be referred as an ’’ideal RG trajectory”: along it one has the exact one-parameter scaling, which is a pipe 
dream in many fields of physics (see e.g. [50]). To define it rigorously, let consider the limit a —>• 0 with 
fixed £/a; then all trajectories lying at the surface £/a = const (Fig.14) tend to one point (analogously to 
the critical surface), while the locus of such points is the ideal RG trajectory. 

Let the parameter k is measuring the distance along the ideal trajectory: then £/a (or A/to) is a 
function of n. Analogously, all dimensionless quantities depend only on k. while the dimensional quantities 
are measured in units of A. As a result, we come to equations 

9 — F g (k) , to = A F m (rc) , = A di Fi (k) , (80) 

which coincide with (76) and give the relations (77) without dependence on to/A. 

The above construction has a following sense. If the limit a — > 0 is taken in the arbitrary manner, then 
the system will go to infinity along the unstable direction and appear far from the critical surface, which 
is our goal. Therefore, we suggest to take the continual limit in two steps: 

(a) take a limit a —y 0 for a/£ = const ; 

(b) take a limit a/£ —>• 0. 

It appears, that dependence on a in relations (77) disappears already at the first step. The second step 
becomes unnecessary and there is no need to take the continuum limit in the bare theory^. The const 

17 These ideas are close to the QCD specialists, and in fact the above consideration was partially taken from ’’Introduction to 
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appearing in (a) is one of the possible definitions of the parameter re. 

Any RG trajectory is a line of ’’constant physics”, since the RG transformation is simply a mental 
construction and does not affect the large-scale properties of the system. All trajectories belonging to the 
critical surface and terminated in the fixed point, give the equivalent lattice representations for the unique 
continuous field theory; another equivalent representation is given by the ideal RG trajectory originating 
from the fixed point. Consider the trajectory AB , which begins near the critical surface and goes along it, 
and then tends to the ideal RG trajectory (Fig.14). Introducing re as a distance along AB , we come to the 
parametric representation analogous to (80) and relations (77), following from it. The latter relations will 
be the same as obtained from (80), since in both cases they are independent of £/a and correspond to the 
physically equivalent models for £/a = 0 and £/a = oo. We can retain definition of re as a distance along 
the ideal trajectory, and assign it to the point of AB, corresponding to the same value of £/a. Writing 
relations analogous to (80) 

g = F g (re) , to = A F m (re) , Ai = A di Fi (re) , (81) 

we see that the second relation remained unchanged, but the rest of them become different. The charge 
g usually belongs to irrelevant parameters and we can introduce ’’the axis of charges” on the critical 
surface; using trajectories of AB type with different directions relative to ’’the axis of charges” we can 
obtain different functions g = F g (n). It means that the functional relation between g and re becomes 
indeterminate and can be omitted. 

As a result, the renormalized and bare sectors of theory become decoupled. The renormalized sector 
contains relations (77), where g and m are considered as independent variables. The bare sector contains 
only relation a/£ = to/A = F m (re), which determines re as a function of a and is irrelevant from viewpoint 
of physics. Parameter a/£ becomes absolutely free. 

Thereby, we come to the following conclusion: renormalizable theory of the considered type allows 
representation in the form of lattice theory, which gives the correct relations between physical quantities, 
and contains free parameter a/£, which does not enter these relations. 

10. Application to theory of confinement 

Consider a pure Yang-Mills theory, i.e. QCD without quarks. Then the quark mass does not enter 
as a parameter and the theory contains no natural mass scale. To avoid the specific difficulties related 
with such situation, we introduce the ’’extended version” of Yang-Mills theory, where the role of the bare 
mass mo (more exactly, the ratio too /A) is played by some auxiliary parameter p characterizing the lattice 
theory; as a renormalized mass, we accept the mass to of the lightest glueball (the bound state of several 
gluons), while the correlation length £ is defined as to" 1 . Thereby, two bare parameters go and p provide 
the observable values for renormalized g and to. In order to return to the standard variant of theory, we 
should remove the introduced extra degree of freedom by fixing one relation between observable quantities. 
However, it can be done on the late stage (see below), while the main of analysis is produced for the 
” extended version”. The latter is analogous to c/> 4 theory. 

According to Wilson [62], confinement can be proved in the lattice version of the Yang-Mills theory for 
large values of the bare charge go. The energy of interaction for two probe quarks separated by a distance 
R is V ( R ) = cr R, while the string tension a and the glueball mass to are given by expressions |611 (521 [55] 

„ = kM, m = ( 82) 

a 1 a 

In spite of the evident success, the Wilson theory is considered as purely illustrative and having no relation 
to real QCD. As was indicated by Wilson himself, his theory corresponds to a situation 

^ Ca or m^> A, (83) 

lattice QCD” by R. Gupta |61| . This picture is discussed there in relation to improvement of the lattice action, and the author 
claims that ’’simulations, done along the ideal RG trajectory, will reproduce the continuum physics without discretization 
errors”. It implies the absence of a/£ dependence, in accordance with our results. Only final conclusion was not made, that 
the continuum limit is not necessary in the bare theory. 
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which is considered as nonphysical. An attempt to advance into the physical region inevitably destroys 
the strong coupling regime. Indeed, the /3-function in terms of the bare charge is believed to be negative 
I64(, 1651 . and go tends to zero in the continuum limit a —> 0. Therefore, the strong coupling regime is 
inevitably destroyed in the physical region £/a > 1 and Wilson’s theory becomes inapplicable. 

A situation changes drastically, if we use representation (76, 77) introduced in the previous sections. In 
this case: 

(i) due to absence of the £/a dependence in (77), this parameter can be taken arbitrary: it eliminates 
objections against the nonphysical regime in Wilson’s theory. 

(ii) there is no direct relation between the bare and renormalized charge; rewriting the second expression 
(82) in the form 



we see that, independently of renormalized values of g and m, it is possible to choose the free parameter 
a/t ; so as to obtain a sufficiently large value for go. Then Wilson’s theory becomes applicable and the first 
relation (82) gives a finite value for a, i.e. confinement. 

Representation (76,77) cannot be introduced for the simplest Wilson action HU EH EH], since it does 
not contain a sufficient number of parameters. To obtain the observable values of a and m one should fix 
both go and a; but the fixed a means impossibility to introduce a representation with free parameter a/£. 
This problem can be easily solved using more complicated forms of the action m- As a result, we obtain 
instead of (82) 

_ _ ln(3gg) _4 ln(3pg) 

cy n • iyb i 

af a 2 

where a± = Aqa, a 2 = k 2 a simply by dimensional reasons [19J . 

Relation (77) in application to a has a form 

a = m 2 F a (g) , 

so g is functionally related with a /to 2 . On the other hand, Eqs. 85 give 

a_ _ a\ 
to 2 16af 1 ii(3<7q) ’ 

and the ratio cr/m 2 changes from a finite value till zero in the strong coupling region go > 1. It means that 
only restricted range of g values can be reproduced. Such restriction is natural due to the physical essence 
of the problem. Indeed, the linear confinement potential is expected only at large distances, where g is 
certainly not small; hence, small values of g are inaccessible in the Wilson regime. Contrary, the restricted 
range of cr/m 2 values goes across the logic of theory. Indeed, a/£ is a free parameter and all physical 
results can be obtained at its arbitrary value. In the case a/£ 1, the regime of confinement is controlled 

analytically and any physically accessible value of cr/m 2 should be possible in this limit. In fact, the range 
of cr/m 2 values can be extended if we use the models with essentially different 01 and 02 m- Absence of 
restrictions on cr/m 2 in the presence of restrictions on g is possible only if dependence cr/m 2 = F a (g) is 
singular (Fig.15); fortunately, we can demonstrate that it is a probable variant. 

Investigations of complicated lattice versions of Yang-Mills theory HI] show existence of phase transi¬ 
tions (lying in the region go ~ 1), corresponding to vanishing of the lightest glueball mass m, with finite 
values of a and other mass parameters. These transitions are considered as lattice artifacts, since they 
do not survive in the continuum limit, when go — ¥ 0. In our approach the limit a — > 0 is not necessary 
and such phase transitions acquire the physical sense. If m vanishes at the point g = g* , then dependence 
alm 2 = F a (g) has a form shown in Fig.15. 

Existence of points with m = 0 in the parametric space means that the ’’extended version” of Yang- 
Mills theory does not possess the mass gap. To eliminate this defect, we should return to the standard 
variant of theory, fixing one relation between the observable quantities. The character of such relations is 


(85) 

( 86 ) 

(87) 
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Figure 15: Dependence of cr/m 2 against g. In order to obtain the special values Cq\ corresponding to 
zero values of the mass gap, one should mark all stable fixed points on the horizontal axis and make a 
construction shown in the figure. 


well-known and is determined by ’’dimensional transmutation” m Sec. 14.1], |6]2 Sec. IV.6], according to 
which all quantities of the same dimensionality differ only by the constant factor, independent of go. For 
our purposes it is convenient to accept the condition 

cr/m 2 = c, (88) 

which defines the one-parameter family of Yang-Mills theories with different values of the structural con¬ 
stant c. Under condition (88), the points with m = 0, a = const become inaccessible. 

It does not yet prove the existence of a mass gap, since er and m can vanish simultaneously. In order to 
analyze such situations, consider the Gell-Mann - Low equation for the renormalized charge g attributed 
to the length scale to -1 

, , ^j = P(g 2 ) = Pog 4 + Pig 6 + • ■ ■ , (89) 

a In md 

where /3-function does not coincide with that of the bare theory, but has the same first coefficients /?o 
and /3i. It is clear that the value g* (Fig.15) is a root of the /3-function; generally, it has several roots 
determining the RG fixed points. In the limit m —> 0, the charge g tends to one of these fixed points, while 
following variants are possible for o/m 2 : (a) cr/m 2 — > oo, (b) a/m 2 —> 0, (c) o/m 2 —>• Co- The first two 
variants are incompatible with Eq. 88, while the third variant is possible in the case c = Co- If there are 
several stable fixed points gP\ then there are several special values (see Fig.15), for which the mass 
gap vanishes; for all other values of c the mass gap is finite. 

Physically, it looks most probable that only one fixed point g* with o /to 2 —> oo is present, so no 
special values arise. Mathematically, one can suggest an infinite number of fixed points, which form 
a sequence Cq 1 '* everywhere dense in the interval (0,oo). However, small values of cr/m 2 correspond to the 
Wilson regime where finiteness of cr and m is verified immediately. As a result, the proof of the mass gap 
is complete for small values of the structural constant c. 

If quarks with a zero mass0 are introduced, then the regime of dimensional transmutation is conserved 

18 In the case of fermions, the renormalization of mass has a multiplicative character and the choice of the zero bare mass 
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and the trick with ’’extension” of theory remains possible; it seems, that the general structure of theory is 
also retained. 


In conclusion, the properties of continuous Yang-Mills theory can be reproduced by a certain lattice 
theory. The bare charge go in this lattice theory can be taken arbitrary, and in particular infinitely large. 
For large go, any reasonable lattice version of Yang-Mills theory gives finite values of er and m. Vanishing 
of m is possible under exceptional conditions, which are avoided in the general situation. As a result, the 
problem of analytical proof of confinement and the mass gap can be considered as solved, at least on the 
physical level of rigor. 


11. Conclusion 

We have given an overview of field theory evolution from its early stage to the present time. A 
widespread opinion on triviality of (f > 4 theory and QED should in fact be revisited. They are trivial in 
Wilson’s sense, but escape the Landau ’’zero-charge” situation. This conclusion is reached by summation 
of weak-coupling perturbation expansions for /3(g), analytical calculation of its strong coupling asymp¬ 
totics, and confirmation of results by analysis of strong coupling expansions. The second possibility in the 
Bogoliubov and Shirkov classification is shown to be actual, allowing to construct the continuous theory 
with finite interaction at large distances. A possibility to retain a lattice in the bare theory allows to justify 
the Wilson approach in theory of confinement. 
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